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iHj ABSTRACT 

This paper presents an analysis of the local peculiar velocity field based on the Wiener Fil- 
^ ter reconstruction method. We used our currently available catalog of distance measurements 

containing 1,797 galaxies within 3000 km/s: Cosmicflows-f . The Wiener Filter method is used 
to recover the full 3D peculiar velocity field from the observed map of radial velocities and to 
recover the underlying linear density field. The velocity field within a data zone of 3000 km/s is 
decomposed into a local component that is generated within the data zone and a tidal one that 
is generated by the mass distribution outside that zone. The tidal component is characterized 
by a coherent flow toward the Norma-Hydra-Centaurus (Great Attractor) region while the local 
component is dominated by a flow toward the Virgo Cluster and away from the Local Void. A 
IL" detailed analysis shows that the local flow is predominantly governed by the Local Void and the 

. !^ Virgo Cluster plays a lesser role. The analysis procedure was tested against a mock catalog, ft 

is demonstrated that the Wiener Filter accurately recovers the input velocity field of the mock 
H catalog on the scale of the extraction of distances and reasonably recovers the velocity field on 

significantly larger scales. The Bayesian Wiener Filter reconstruction is carried out within the 
ACDM WMAP5 framework. 

The Wiener Filter reconstruction draws particular attention to the importance of voids in 
proximity to our neighborhood. The prominent structure of the Local Supercluster is wrapped 
in a horseshoe collar of under density with the Local Void as a major component. 

Subject headings: galaxies: distances and redshifts, cosmology: observations- large-scale structure of 
Universe 
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1. Introduction 

This paper presents a reconstruction of the 
large scale structure of the nearby Universe from 
a survey of galaxy peculiar velocities using the 
means of the Wiener Filter (WF) and constrained 
realizations (CR). The observational component 
is a catalog of distances and inferred line-of-sight 
peculiar velocities: Cosmicflows-1, data assembled 
by Tully et al. (2008). This information can be 



used to derive three-dimensional velocity and den- 
sity fluctuation fields with use of the Wiener filter 

The subse- 



Zaroubi et al. 


1995 


Hoffman] 


2009 



quent description of phase space in the local region 
can be used to infer initial conditions that can be 
the starting point for constrained numerical simu- 
lations dGottlober et al.||2010[ ). 

The goal of the program is to strengthen the 
link between an improving knowledge of the ob- 
served positions and motions of nearby galaxies 
and the increasingly sophisticated theoretical ren- 
ditions of the development of large scale structure 
and galaxy formation. It is exquisitely well es- 
tablished from the observation of the dipole tem- 
perature variation of the Cosmic Microwave Back- 
ground that our Galaxy has a deviant motion of 



630 km/s ( Fixsen et al.||1996 ) but there is a very 
poor understanding of how that arises. There are 
local components like our collision course toward 
M31( Hoffman et al. 2007), the repulsion by the 



Local Void, and the pull from the Virgo Cluster 
and surrounding structure ( ,Tully et al.||2008 ). On 
a larger scale, it is now apparent that the Galaxy 
lies at the periphery of an overdense region cen- 
tered near the giant clusters in Norma, Hydra, and 
Centaurus at distances of 3000—5000 km/s. It has 
been long established that we are part of a flow in 
that direction ( Lilje et al.|1986 Lynden-Bell et al 
1988). However there has been a lively debate 
about whether the flow ends there, or continues 
on to the Shapley Concentration at 15,000 km/s 



(Pike & Hudson 


2005 


Kocevski & Ebeling 


2006 


Erdogdu et al. 


2006 


Lavaux et al. 


20101 ] 


Bilicki 


et al.,|2011) , or continues even further (Feldman 


et al. 2010 Kashlinsky et al. 2010 


. There are 



important things to be learned from cosmic flows. 

The discussion will begin with a description of 
current observational material, with brief men- 
tion of anticipated developments. This section 
will be followed by a discussion of Wiener Filter 



reconstruction of three-dimensional velocities and 
a density field compatible with the observations 
within the limitations of a linear analysis. The 
viability of the Wiener Filter constrained realiza- 
tions is checked in the next section with mock cat- 
alog comparisons. A constrained simulation that 
provides a model that is in general but not fine- 
tuned accordance with the local structure provides 
the test. It will be demonstrated that the proce- 
dure does a good job of recovering the observed 
structure and velocity field in the domain of obser- 
vations and does a reasonable job of anticipating 
the structure at distances well beyond the obser- 
vations. 

2. Data 

This article mixes observations and theoreti- 
cal modeling. Observed distances are based on 



a generally accepted zero-point scale (Freedman 



et al. 2001 ) so are reported in Mpc. The nat- 
ural distance unit for the models is Mpc 
where h expresses the Hubble Constant Hq = lOO/i 
km/s/Mpc. Distances related to interpretation of 
models are often quoted as h^^ Mpc and h = 0.74 
gives an appropriate translation to the observa- 
tions. If a distance is not tagged with the param- 
eter 'ft,' then it is on the observed scale. 

2.1. Observed peculiar velocities 

In |Tully et al] ( |2008[ ) we published a cata- 
log of distances that will be hereafter identified 
as "Cosmicflows-l'l^ The catalog contains dis- 
tances for 1797 galaxies within 3300 km/s, draw- 
ing from four different methods: the Cepheid 
Period-Luminosity relation. Tip of the Red Gi- 
ant Branch and Surface Brightness Fluctuation 
luminosity indicators, and the galaxy luminosity- 
rotation rate or TuUy-Fisher correlation. The first 
three methods give distances with accuracies of 
10%, but in the case of the first two only to very 
local galaxies, and in the case of the third only 
to early types that are preferentially located in 
high density regions. The Tully-Fisher method, 
invoking the correlation between the neutral HI 
gas rotation rate and the integrated magnitude of 
a spiral galaxy, has an rms accuracy of 20% with 



^Available in the Extragalactic Distance Database I Tully 
|et al.|2009| ^http:/ /edd.ifa.hawaii.edu^ 
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measurement in a single band but the lesser accu- 
racy per measurement is offset by the availability 
of large samples widely distributed over the sky 
and in distance. 

Peculiar velocities {Vpec in km/s) and corre- 
sponding errors are computed from distances (d 
in Mpc) and their errors assuming a value for the 
Hubble constant Hq = 74 km/s/Mpc ( |Tully et al. 



2008) 



pec 



CMB 



d * Hq (km/s) 



-Ad*Ho (km/s) 



(1) 
(2) 



The distribution of rms percent errors in dis- 
tance is shown in Figure [ij The peak at 20% re- 
sults from single measures with the Tully-Fisher 
method. Distance measures with uncertainties 
larger than 20% are not accepted as sufficiently 
accurate and thus are not used in this catalog. As 
seen on Figure [T] the median rms error in both 
distance and pecuhar velocity is 13%. 

It is to be appreciated that the Cosmicflows-1 
catalog is not a complete sample. Incompleteness 
does not affect the reconstruction of the three- 
dimensional velocity field as long as the underly- 
ing individual distances are not subject to biases in 
the distance measures]^ The observed peculiar ve- 
locity field can be sparse in an inhomogeneous way 
without leading to deviant results in the recon- 
struction as it will be demonstrated later in this 
paper. Figure[2] shows the distribution of galaxies 
with distance measures in the Cosmicflows-1 cat- 
alog in supergalactic coordinates SGL, SGB. The 
velocity distribution of the sample is shown in Fig- 
ure [T] 

The general features of the large scale structure 
within 8000 km/s are illustrated in Figure |2] This 
figure is a representation of the 30,124 galaxies in 
the redshift catalog drawn from the literature that 
we call 'V8Kj^ The 1,797 data points with mea- 
sured peculiar velocities within 3,000 km/s are in- 
cluded in this larger catalog. The WF is sensitive 
to structure on scales larger than the domain of 
the distance measurements. The V8K catalog is 



^The subject of biases is discussed in detail by (TuUy &; 
|Pierce|2000| |. 

*The V8K catalog is available at http://edd.ifa.hawaii.edu 



used in this paper to compare the observed and 
derived densities from small to large scale^ 

The catalog of distances is a work in progress. 
Data relevant to the measurement of distances are 
being accumulated in the Extragalactic Distance 
Database, and an extended version of this catalog 
" Cosmicflows-2" is in preparation. Observations 
are being extended to ~ 15,000 km/s in a pro- 
gram of HI observations and photometry that we 
are calling " Cosmic Flows" ( Courtois et al.|20TT I . 
HI observations are being made in the context of 
the Cosmic Flows Large Program with the 100m 
Green Bank Telescope at the National Radio As- 
tronomy Observatory and complementary south- 
ern observations with the Parkes Telescope in Aus- 
tralia. HI profiles for more than 13,000 galax- 
ies have been measured consistently by the same 



pipeline (Courtois et al. 20091 ). The photometry 
for the program is being acquired from the Uni- 



versity of Hawaii 2.2m Telescope (Courtois et al 



2011), the multiband survey Pan-STARRS in the 



North and, it is anticipated, SKYMAPPER in the 
South. 

3. Analysis: Theoretical Tools 

The analysis of the velocity field is performed 
here within the standard model of cosmology. 
The model assumes that structure has evolved 
out of a primordial Gaussian random perturbation 
field via gravitational instability. This hypothe- 
sis implies that on large enough scale the flow is 
curl-free and that the velocity and gravitational 
fields are linearly related. The Bayesian Wiener 
Filter Constrained Realization (WF/CR) analy- 
sis is performed here within the standard ACDM 
model assuming the cosmological parameters de- 



rived from the WMAP5 data base (Komatsu et al 



2009 ) . The method is tested against mock catalogs 



drawn from a constrained simulation run within 
the ACDM/WMAP3 framework. 

3.1. Wiener Filter and Constrained Real- 
izations 

Peculiar velocity surveys constitute a very noisy 
and sparse data set diluted across the whole sky. 



^Alternatively, comparisons could be made with 2MRS, the 
2MASS Redshift Survey (Huchra ct al. 2005) which extends 
farther. The two catalogs arc qualitatively similar but V8K 
includes fainter galaxies so provides denser coverage. 
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Nonetheless, the pecuhar velocities directly sam- 
ple the gravitational field induced by the large 
scale structure (LSS) of the universe and there- 
fore can be used to uncover the underlying den- 
sity and velocity field. The WF algorithm pro- 
vides the optimal tool for the reconstruction of 
the LSS from peculiar velocity surveys. The gen- 
eral WF framework has been thoroughly reviewed 
inlZaroubi et al.|(|1995|) andlHoffmanl (|2009|). The 



application of the WF to peculiar velocity obser- 



vations was first presented in Zaroubi et al. ( 1999 1 



In the standard model of cosmology the LSS 
emerges out of a primordial random Gaussian per- 
turbation field via gravitational instabilities. In 
the linear regime of deviations from the homo- 
geneous and isotropic Friedmann expansion the 
density and velocity fields retain their Gaussian 
nature. The WF constitutes the optimal linear 
minimal variance estimator given a general data 
set and an assumed prior cosmological model. It 
can be shown that in the case of an underlying 
Gaussian random field the WF coincides with the 
optimal Bayesian estimator of the field given the 



data (Zaroubi et al. 1995) 



The linearity assumption implies that the den- 
sity and velocity fields are reconstructed as if the 
linear regime is valid today on all scales. How- 
ever, the present epoch LSS is non-linear, with a 
smooth transition from the very large scales where 
the linear theory prevails down to small scales that 
are governed by non-linear dynamics. The scale of 
transition from the linear to the non-linear regime 
depends on the nature of the field under consid- 
eration. The scale of non-linearity of the velocity 
field is roughly a few Mpc and for the density field 
it is about a factor of ten larger. It follows that 
the WF algorithm provides an excellent tool for 
the reconstruction of the present epoch velocity 
field. 

A final comment on a property of the WF recon- 
struction is due here. The WF constitutes a con- 
servative estimator which is designed to find the 
optimal balance between the observational data 
and the prior model one is assuming in construct- 
ing the WF (Zaroubi et al. 1995). In regimes 



where the data is dense and accurate, the WF 
is dominated by the data while in regimes where 
the data is very noisy or very sparse the WF re- 
covers the prediction of the prior model. In the 
standard model of cosmology the expected mean 



density departures and velocity fields are zero. In 
observational surveys of peculiar velocities the sta- 
tistical error of individual measurements increases 
linearly with distance and the data sparseness in- 
creases with the distance. The degradation in the 
quality of the data with the distance implies that 
the WF reconstruction is expected to attenuate 
the reconstructed fields towards the null fields. 

The WF presents the conditional mean field 
given the data and the assumed prior model. The 
variation of the actual field around the WF mean 
field can be sampled by performing CRs of the 



underlying Gaussian density field ( Bertschinger 



1987 Hoffman & Ribak 1991 Zaroubi et al. 1999 1 



The scatter of the CRs provides a measure of the 
statistical uncertainties of the WF reconstruction. 

3.2. Decomposition into Tidal and Local 
Flow Fields 

An interesting question given attention by mod- 
elers of the velocity field concerns the identifica- 
tion of the sources causing the motion. What cos- 
mological objects such as clusters, superclusters 
and voids, are affecting the velocity field within a 
given region of space? 

For a given reconstruction of the 3D velocity 
field in a given region of radius R, the velocity 
field can be decomposed into two components, one 
that is induced by the mass distribution within 
the volume and one induced by the outside mass. 
This decomposition can be easily done in the lin- 



ear regime (Hoffman et al. 2001) where the ve- 
locity field is proportional to the density field and 
therefore the velocity and density fields are related 
by a Poisson-like equation: 



(3) 



Given a density field within a given volume and as- 
suming boundary conditions on that volume, the 
Poisson-like equation can be solved and the veloc- 
ity field can be calculated. This amounts to the 
particular solution, to which one can always add 
the homogenous solution of the source-free equa- 
tion, namely the Laplace equation. The conse- 
quence is the separation of the velocity field inside 
a given volume into its local component, namely 
the one induced by the internal mass distribution. 



and the tidal component (Hoffman et al. 20011 



It is to be noted that we used the term 'tidal' 
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to represent the full velocity field that is induced 
by the external mass distribution, including the 
dipole term (i.e. the bulk motion of the volume 
due to external causes at any depth). 

The WF reconstructs the density and velocity 
fields on a grid that spans some computational 
box. Alternatively, one can use the WF veloc- 
ity field and calculate the linear density field by 
taking the divergence of the velocity field. Sup- 
pose that the velocity field is to be separated into 
the component induced by the density within a 
sphere of radius R, enclosed within the box, and 
the residual of that component. This separation is 
done by solving Eq. |3]for a density field set to the 
WF field within R and taken to be zero outside R. 
In practice the operation is carried out using an 
FFT solver and imposing zero-padding outside the 
desired volume. The resulting field is the diver- 
gent or local component, and the residual of the 
full WF velocity field after vector subtraction of 
the local component constitutes the tidal field in- 
duced by the mass distribution beyond the sphere 
of radius R. Visual inspection of the local field 
often suffices to infer the existence of "great at- 
tractors" that lie outside the sphere. A multipole 
decomposition of the tidal field provides a more 
quantitative characterization of the tidal field. 

3.3. Structure Outside the Data Zone 

Peculiar velocity data sets consist of a finite 
number of data points taken within a finite vol- 
ume, the data zone. Peculiar velocities have a 
much larger correlation length than the underly- 
ing density field, ft follows from the latter fact 
that the WF reconstruction can be applied on a 
volume larger than the data zone. The extent and 
quality of the reconstruction depends on the as- 
sumed prior model, i.e. the power spectrum, and 
the quality of the data: the magnitude of the er- 
rors and the sparseness and the space coverage of 
the data. In the following, the WF reconstruc- 
tion is to be conducted within two boxes of size 
i = 64 and 160h^^Mpc, centered on the origin of 
the Supergalactic coordinates system. 

3.4. Testing with Mock Catalogs 

A detailed analysis of the tests conducted with 
mock catalogs is available in the accompanying pa- 
per (Hoffman et al. 2011). We give an overview 



here. A mock velocity catalog has been drawn 
from a constrained simulation of the local uni- 
verse. The simulation was performed within the 
CL UES projeclpj specifically the case discussed 
by Cuesta et al. (2011 1. The simulation is con- 
tained within a box of L = 160h~^Mpc, hereafter 
BOX160, assuming the WMAP3 cosmological pa- 
rameters. This model is a revision of the |Klypin| 



et al. (2003) constrained simulation, with a differ- 



ent realization of the random field and more mod- 
ern estimation of the cosmological parameters. A 
visualization of the simulation and some discussion 
of the structural properties is presented in |Cuesta| 



et al. (2011 1. The simulation reproduces the main 
players of the local cosmography. In particular 
it recovers a representation of the Local Super- 
cluster, the Virgo and Coma clusters, the Great 
Attractor and the Perseus-Pisces Supercluster. It 
also reproduces an environment like the neighbor- 
hood of the Local Group (LG) and its cold Hubble 
fiow. The simulation provides a very good dynam- 
ical replication of the local universe and thereby 
constitutes a good 'laboratory' for tests of our 
tools of analysis. 

This is a dark matter only simulation. A hier- 
archical friends-of-friends halo finder was applied 
to the simulation and a halo catalog has been con- 
structed. A mock galaxy catalog has been con- 
structed using the algorithm of the conditional 
luminosity function (van den Bosch et al. 20071. 



An 'observer' has been placed at the position of 
the LG and the peculiar velocity of 1,797 galaxies 
has been measured, assuming a Poisson distribu- 
tion of errors centered around 13%. The selection 
of the galaxies is done randomly in a sphere of 
R — SOh^^Mpc centered on the LG in order to 
mimic the observed redshift distribution. 

The upper right panel of Figure [3] shows the 
Supergalactic Plane of the constrained BOX160 
simulation. The density and velocity fields are 
Gaussian smoothed with Rg = 1 h~^Mpc. At such 
high resolution the density field lies well within 
the non-linear regime. The WF reconstruction is 
applied to the simulation with alternative smooth- 
ing lengths oi Rg — 5 and 10h~^Mpc. A visual 
inspection of the maps reveals that the reconstruc- 
tion recovers the overall LSS of the box. The com- 
parison in Table [T] shows that the reconstruction 

^http:/ /www. CLUES-project.org 
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V{Rg = 5h-^Mpc) 


V{Rg = lOh-^Mpc) 


Slope 


0.81 


0.92 


Corr. Coeff. 


0.89 


0.94 


Scatter 


57 km/s 


35 km/s 


Mean 


9.2 km/s 


3.5 km/s 


St. Dev. 


63 km/s 


36 km/s 


Skewness 


0.41 


0.26 


Kurtosis 


0.85 


1.12 



Table 1: Comparison of velocity fields with the WF reconstructed fields from the mock catalog. The velocity 
field is represented by the magnitude of the velocity vector. The comparison is done for Rg = 5 and 
10h~^Mpc. The first three rows present the slope, correlation coefficient and scatter of the linear regression 
of the simulation and the WF reconstruction. The following rows give the first 4 moments (mean, standard 
deviation, skewness and kurtosis) of the residual of the WF reconstructed fields from the simulated fields. 



is best with the larger Gaussian smoothing and 
is excellent within the 30h~^Mpc annulus of the 
data zone. The extrapolation at larger distances 
recovers the gross features of the LSS but at a 
degraded quality. It follows that the WF can be 
used to provide a rough expectation of the LSS in 
unobserved regions. 

A linear regression of the N-body and WF re- 
constructed fields is given in Table [l] The table 
also provides the mean, standard deviation, skew- 
ness and kurtosis of the residual of the WF from 
the N-body fields. 

4. Cosmography 

4.1. The volume within 40 Mpc 

The series of panels within Figure |4] illus- 
trate the broad features of the WF reconstruction 
within the radius containing distance measure- 
ment constraints; the data-zone. The contours of 
colors illustrate the WF reconstructions based on 
the observed velocity field. The left panels display 
linear density and the right panels display radial 
peculiar velocities. In the left panels the individ- 
ual points that are plotted are galaxies drawn from 
the V8K catalog that lie within the specific dis- 
tance shells. The pronounced correlation of galax- 
ies with the high density WF contours provides 
a visual demonstration that the WF methodol- 
ogy can recover the 3D density distribution of the 
local universe from a sole knowledge of peculiar 
radial velocities. It is to be noted here that the 



V8K galaxies are a rough proxy for the actual 
spatial distribution of nearby galaxies. Galaxies 
are plotted based on their redshifts, assumed to 
be rough estimators of their distances. The WF 
reconstructed radial velocity field is shown in the 
right panels. 

Figure [5] provides orthogonal views of slices of 
the WF reconstruction. The panels on the left il- 
lustrate densities as contours and velocities as vec- 
tors. Model velocities are black and on a grid and 
observed radial velocities are in color (red: out- 
ward; blue: inward). The panels on the right show 
densities again but now with color and galaxies 
within the slices from the V8K redshift catalog are 
superimposed. With the top two rows the slices 
in the normal direction are on the cardinal axes 
(SGZ==0 and SGX=0 respectively) while with the 
bottom plots the shce is displaced to SGY=16.9 
Mpc to coincide with the main structure of the 
Local Supercluster which runs through the Virgo 
Gluster on a line roughly parallel to the SGX axis. 
This structure is seen top-down in the upper pan- 
els and edge-on in the lower panels. 

The WF reconstruction only provides a descrip- 
tion of the linear regime and hence it has a lim- 
ited power to reconstruct the present epoch struc- 
ture on small scales. However, the WF is shown 
to be very good for a recovery of the important 
large scale features. Perhaps the most interest- 
ing revelation with the current study is the nature 
of the under dense region that wraps around the 
Local Supercluster centered on the Virgo Gluster 
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(SGY=17 Mpc near SGX = SGZ = 0). The larger 
part at SGZ posit ive and -20 < SGY < +10 Mpc 
is the Local Void ( jTully fc Fisher|1987|[Tirily et al 
2008 ) . There is also a pronounce void directly be- 



hind the Virgo Cluster at SGY - +30 Mpc and 
SGX ~ that we will call the 'Virgo Void'. It be- 
comes evident here (best seen in the middle panels 
of Fig. [5]) that the Local and Virgo voids are linked 
by a generally under dense region that raps above 
the Local Supercluster at SGZ ^ 15 Mpc. Sparse 



filaments do course through this volume (TuUy & 



Fisher 1987). However the overriding aspects of 
our local region are, first of all, the general under 
density across a region of ~ 60 Mpc but, secondly, 
the bridge of over density of the Local Superclus- 
ter at SGY ~ 18 Mpc, SGZ ~ and extending to 
negative SGZ, and a width in SGX that establishes 
links to major structures including the Centaurus 
and Hydra clusters (seen best in the top panels of 
Fig.[5]). 

These features are seen in other reconstruc- 



tions of the local neighborhood ( 


di Nella et al. 


1996 


Branchini et al.||1999 


Erdogdu et al.||2006|. 


Romano-Diaz & van de Weygaert 


(2007 


) empha- 



et al. (2009) provide hints of voids on scales of 
150 Mpc. The present study is distinguished by a 
higher density of distance information on a smaller 
scale and, in particular, benefits from the observa- 
tion of the kinematic signature of void expansion; 
the discontinuity in peculiar velocities between 
galaxies in the local wall of the void and galaxies 
beyond the wall ( |Tully et al.|2008[ ). The s ame dis- 
tance constraints were used in the study by |Lavaux| 
et al. (20101 but only to confirm the properties of 



the bulk flow of the local region within 3000 km/s. 

The series of panels in Figure [6] present the evo- 
lution of the local field on Aitoff all-sky projections 
at increasing distances. The plotted velocity field 
is the radial part only of the local field associated 
with mass within 40 Mpc. The top panel at 6.7 
Mpc shows a motion from Supergalactic north to 
south and toward the direction of the Virgo Clus- 
ter, which is beyond this slice. The panels at 13.5 
and 20.3 Mpc display the front fall away from us 
toward the Virgo Cluster and backside infall to- 
ward us. The same effect is seen around the For- 
nax Cluster: a front and a backside infall around 
sgl=270, sgb=— 45. The color bar gives the veloc- 
ities in km/s. Within this local volume the largest 



motion is toward a location at sgl=150, sgb=0. 
The farthest shells show the fading of the WF re- 
construction towards the null fields as it should 
from construction. 

4.2. Extrapolation outside the data-zone: 
the Great Attractor 

The WF reconstruction has two modes of oper- 
ation. On the one hand it interpolates the struc- 
ture in between the data points, providing the op- 
timal interpolation given the data, its errors and 
the underlying prior model. On the other hand 
the WF extrapolates the structure to regions not 



covered by the data ( 


Zaroubi||1995 


Zaroubi et al. 


1995 


Hoffman et al. 


2001 


I. The long range co- 



herent nature of the peculiar velocities makes the 
WF an effective tool to probe large scale struc- 
ture beyond the spatial confines of the data. In 
the current situation, the catalog of observed pe- 
culiar velocities extends to 3,000 km/s. In order to 
probe to larger scales the WF has been applied to 
a 160h~^Mpc box, centered on the Local Group, 
on a 128"^ grid. 

Results on this larger scale are illustrated in 
Figure [7] The top left and right panels are ex- 
panded scale versions of the top and middle pan- 
els of Fig. [5] As was described in Section 3, the 
WF density reconstruction degrades to the mean 
density on scales beyond the data-zone. The lo- 
cal component of the reconstructed velocity field is 
dominated by the repulsion of the voids and fiows 
toward the densest structure. The dominant con- 
vergence at SGX 40 Mpc, SGY ~ +20 Mpc 

is coincident with prominent structure in Norma, 
Hydra, Centaurus, and Pavo-Indus, the principal 
components of what came to be called the 'Great 
Attractor' ( |Shaya|1984 [ [ Dressier et al.|1987[ ). Con- 
sideration of the tidal part gives the visual impres- 
sion that the large scale flow on a scale of 80 Mpc 
is dominated by attractions outside the computa- 
tion box. 

Alternative views are given in Figure [8| The 
top panels are the same as with Fig. [7] minus the 
superposition of the V8K galaxies. The lower pan- 
els show the gridded local and tidal decomposi- 
tions respectively, with a decomposition radius of 
80 Mpc. The choice of a decomposition radius is 
very large. It comfortably includes the canonical 
Great Attractor region, located in the region of 
velocity convergences in the middle panels. Even 
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V [km/s] 


Igal 


bgal 


sgl 


sgb 


^bulk 


378 


293° 


18° 


156° 


-21° 


Ai X i? 


261 


233° 


2° 


124° 


-80° 


A2 X i? 


-26 


144° 


0° 


1° 


-7° 


A3 X i? 


-235 


231° 


-88° 


270° 


-8° 



Table 2: Fitting the WF tidal field by a bulk and shear model as in Equation |4j The ampHtude, Galactic 
{I, b) and Supergalactic (sgl, sgb) direction of the bulk velocity and the eigenvalues and eigenvectors of the 
shear tensor are presented. The eigenvalues are multiplied by i? = 30h~^Mpc. 



with this large decomposition radius the tidal flow 
shows a convergence of the flow just beyond this 
radius, thereby indicating the existence of struc- 
ture beyond the Great Attractor. It also suggests 
that the reconstructed Great Attractor backside 
infall arises from the lack of data on this scale and 
is not physical. 

The hypothesis that the tidal velocity field is 
predominantly induced by a Great Attractor con- 
vergence can be tested. The tidal velocity field 
has been fitted by the bulk and shear fiow model 
of Hoffman et al. (2001). Consider Eq. 6 from 



that paper: 



Vair) = Ba + (iiSaP + 'Safj^p 



(4) 



where Be, is the bulk velocity, Hap is the traceless 
shear tensor, HSa/s is a possible local isotropic per- 
turbation, and the indices a and /3 range over the 
3 Cartesian coordinates. The bulk and shear cor- 
respond to the dipole and quadruple terms of po- 
tential theory. The isotropic term H5afj vanishes 
in an analysis of tides. 

Table [2] presents the bulk velocity, Vbuik, the 
three eigenvalues of the shear tensor, Ai, A2 and 
A3, and the direction of their corresponding eigen- 
vectors. The eigenvalues are multiplied by _R = 
30h~^Mpc, so they acquire the dimensionality of 
velocity. A tidal flow induced by a single attractor 
is expected to have: 

a) the eigenvector corresponding to Ai, hereafter 
the first eigenvector, and the bulk velocity to be 
aligned; 

b) the shear eigenvalues have the following struc- 
ture: 

A2 = A3 = — Ai/2. 

The results of Table [2] imply that the local tidal 
flow cannot be modeled by a single attractor. The 



61° misalignment between the flrst eigenvector 
and the bulk velocity is one indicator. Yet, the 
strongest evidence comes from the structure of the 
eigenvalues, which are closer to a (A,0, — A) struc- 
ture which is expected for a random field point 
( [Hoffman et"aI]p0T| ). 



4.3. Bulk velocity 

The bulk velocity of any velocity field is defined 
here as the volume weighted average velocity of 
the full 3D velocity field in a given volume. For 
this discussion, a sphere of radius R — 30h^^Mpc 
is assumed and the velocity field is assumed to 
be evaluated on a regular Cartesian grid. Table [3] 
presents the bulk velocity, Vjjuik, of the the WF re- 
constructed velocity field. To study the robustness 
of WF estimator of Vbuik a suite of 16 CRs have 
been constructed, using the same observational 
data and prior model as in the WF reconstruc- 
tion. The construction of a CR involves the com- 
plementary use of an unconstrained realization of 
the underlying field. These random realizations 
(dubbed RANs) are used as a tool for calculat- 
ing the mean and scatter of the amplitude of the 
bulk velocity of a random sphere in the assumed 
ACDM (WMAP5) cosmology. Both the CRs, and 
therefore also the RANs, are evaluated on a box 
of L = 640h-iMpc on a 128^ grid. This large 
computational box is used to avoid suppression 
of the residual power contributed by very large 
scales. The definition of Vi^uik used here differs 
from the usual definition adopted in the literature. 
Whereas Vbuik is often taken to be the mean veloc- 
ity of a given sample of galaxies, the bulk velocity 
is defined here by the mean velocity of all galaxies 
within a given volume. 

Table |4] provides information on the sources of 
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^huujkm/^] l,,al f>yal 



WF 415 295° 19° 

CRs 401 ± 19 299° ±6 24° ± 4 
RANs 341 ± 130 



Table 3: Statistics of the bulk velocity of the full velocity field. Row 1: magnitude and direction of Vbuik 
in the WF reconstructed velocity field. Row 2: The mean and standard deviation of Vbuik taken over an 
ensemble of 16 CRs. Row 3: The mean and standard deviation of Vbuik taken over 16 random, unconstrained 
realizations. Note that for the random realizations the direction of the bulk velocity is random, hence its 
mean and scatter are not shown. 



V Igai bgai sgl sgb Vy alignment Vcmb 

[km/s] [deg] [deg] [deg] [deg] [km/s] [km/s] [km/s] [deg] 



CMB 


627 


268 


27 


138 


-38 


-364 


330 


-390 





Vlg 


593 


273 


29 


139 


-34 


-374 


326 


-324 


5 


TIDAL 


382 


298 


16 


161 


-17 


-348 


118 


-106 


30 


LOCAL 


304 


234 


36 


98 


-47 


-28 


208 


-220 


30 


VIRGO (7 Mpc//i) 


29 


277 


75 


101 


-3 


-5 


28 


-1 


49 


VIRGO (10.6 Mpc/h) 


60 


202 


18 


45 


-55 


24 


24 


-49 


2C 
60 



65 (w.r.t. Virgo) 



Table 4: Dissecting the Local Group motion. CMB: The CMB dipolo used here to transform the data to the 
CMB frame. Vlg- The LG motion reconstructed by the WF. TIDAL: The Vlg component induced by the 
tidal field, namely the mass distribution outside a sphere of R = 30 h~^Mpc. LOCAL: The Vlg component 
induced by the local mass distribution, i.e. within R = 30h~-'^Mpc (this is the divergent component). VIRGO 
R = 7h~^Mpc: The Vlg component induced by the mass distribution in and around the Virgo cluster, 
specifically within a sphere of radius R = 7h-iMpc centered on [SGX, SGY, SGZ] = [3, 10, 0] h-^Mpc. 
VIRGO R = 10.6h-iMpc: Same as Virgo R = 7h-iMpc but for a sphere of radius R = 10.6h-iMpc, 
i.e. the LG- Virgo separation. Columns 2-4: the absolute velocity and its Galactic (1, b) direction; columns 
5-9: the Supergalactic direction and decomposition of the velocity; column 10: the alignment with the CMB 
dipole. 
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the motion of the LG with respect to the CMB 
frame of reference. The full reconstructed velocity 
field yields a velocity Vlg that is very close to the 
CMB dipole motion, Vcmb, with only 5° misalign- 
ment. The tidal and local components have sim- 
ilar amplitudes and deflections but are mutually 
misaligned by 60° . The last two rows of Table |4] 
present the values of the velocity field at the po- 
sition of the LG induced by the mass distribution 
in spheres of radius 7 and 10.6 Mpc/h centered 
on the Virgo Cluster. The mass in the smaller 
sphere induces a motion closely aligned with the 
cluster but the amplitude of 29 km/s is small. The 
mass in the larger sphere induces a velocity with 
an increased amplitude of 60 km/s but the direc- 
tion is 65° away from the cluster. It follows that 
the Virgo Cluster plays a minor role in shaping 
the motion of the Local Group with respect to the 
CMB. The growth of the amplitude and the shift 
in direction with the increase in radius of the shell 
around the cluster are consistent with the propo- 
sition of a dominant role of the Local Void, as 
discussed in the next section. 

5. Discussion and Summary 

The WF methodology has been applied to the 
Cosmicflows-1 data base which consists of 1,797 
galaxies within 3000 km/s with distance measure- 
ments to reconstruct the large scale structure of 
the local universe. The full linear density and 
3D velocity fields have been reconstructed, assum- 
ing as a prior model the standard ACDM model 
with the cosmological parameters derived from the 
WMAP5 data base. The simultaneous reconstruc- 
tion of both the density and velocity fields has en- 
abled a detailed and critical examination of the 
sources that induce the observed velocity field on 
scales of a few tens of Megaparsecs. The main 
progress since the first application of the WF/CR 
application to velocity data by Zaroubi, Hoffman 
and Dekel (1999) has been along three main lines: 
better observational data, improved testing by 
mock galaxy catalogs, and the employment of a 
tidal field decomposition. A brief summary of our 
main results follows. 

Maybe our most striking finding is the domi- 
nant role of the Local Void in shaping the local 
velocity field. Voids are a prominent feature of 
our neighborhood and the Local Void is of spe- 



cial note. TuUy et al. (2008) concluded that this 



feature extends to SGZ > 40 Mpc and expansion 
of the wall at our location is a significant source 
of the peculiar velocity of our Galaxy. Indeed the 
Local Void is a prominent feature of the WF re- 
construction. The motion identified as 'Local' in 
Table |4] is in excellent agreement in both ampli- 
tude and direction with the local motion identi- 
fied by Tully et al. (a motion of the Local Sheet of 
323±25 km/s toward Igai = 220±7, bgai = +32±6 
or sgl = 80, sgb = —52). The dominance of 
the Local Void marginalizes the dynamical role of 
the Virgo cluster on the local scene. Indeed, the 
structure within a Virgocentric sphere of a radius 
of i? = 10.6h~^Mpc contributes a mere 60 km/s 
which is directed 65° away from the direction of 
the Virgo cluster. It follows that the Virgo cluster 
does not dominate even that Virgocentric sphere. 
Technically speaking, the current analysis is done 
within the linear regime and the results can be 
strictly implied only with that framework. How- 
ever, the dynamical role of a virialized cluster at 
a distance of > 10h~^Mpc away is only modestly 
affected by non-linear effects. 

The motion of the LG with respect to the 
CMB is recovered by the present reconstruction to 
within 5% in amplitude and 5 degrees in direction. 
Decomposing the velocity field into local and tidal 
components, with respect to the data zone sphere 
of i? = 3000h~^Mpc, we find that the local struc- 
ture contributes a component of 304 km/s leaving 
a residual tidal component of 382 km/s, both de- 
viating from the CMB dipole by 30° but in almost 
orthogonal directions. 

Next, the bulk velocity of the full data-zone 
has been examined. The bulk velocity of the 
full velocity field is very robustly determined to 
be Vbuik = 401 ± 19 km/s in the direction of 
Igai = 299 ± 6° and bgai = 24 ± 6°. The sta- 
tistical scatter is calculated over an ensemble of 
CRs. In principle, the bulk velocity can be par- 
tially induced by internal sources although the ex- 
pectation is that it should be dominated by the 
exterior mass distribution. Indeed, the bulk ve- 
locity of the tidal field is found to be 378 km/s 
in the direction of Igai ~ 293° and bgai = 18°. 
The shear tensor component of the velocity field 
does not have the characteristics that would be in- 
duced by a single attractor and none of the eigen- 
vectors are closely aligned with the bulk velocity. 
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It should be noted that the bulk velocity is defined 
here as the volume averaged velocity field within 
the region under consideration. This differs from 
the commonly defined bulk velocity, which is taken 
to be some weighted mean velocity of the galax- 
ies that constitute a peculiar velocities data set. 
It follows that a direct comparison with the bulk 
velocity determinations of other surveys and by 
different methods should be done with care. Hav- 
ing stating this caveat we compare our estimated 
bulk velocity of the R = 30h~^Mpc sphere with 



Nusser & Davis (2011) who studied the SFI- 



catalog and estimated the bulk velocity of a sphere 
of R = 40h-iMpc to be 333 ± 38 km/s toward 
Ig^i = 276 ± 3° and bgai = 14 ± 3°. Given the very 
different methods of estimation and the different 
depth of the studies the two determination are in 
a good qualitative agreement. 

The WF/CR algorithm presented here provides 
a robust methodology for recovering the underly- 
ing density and 3D velocity fields from observa- 
tional surveys of the peculiar velocities of galaxies. 
Further progress in this field is to be pursued along 
two parallel tracks. On the one hand the method 
will be applied to a much larger and deeper sam- 
ple of peculiar velocities measured in our 'Cosmic 
Flows' program ( [Courtois et al.|2011[ ). This effort 
is to be matched by further theoretical work on 
the improvement of the method. The main short- 
coming of the WF/CR method is the fact that 
it is strictly formalized within the linear regime. 
Our work-in-progress focuses on two main issues. 
One is the quest for a good grouping algorithm of 
the data as a mean for suppressing virial velocities 
in collapsed objects and thereby 'linearizing' the 
data. The other is a dynamical mapping that en- 
ables the reconstruction of the quasi-linear density 
field from the WF linearly estimated field. 
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Fig. 1. — Observational data in Cosmicflows-l catalog. Top left. Distance distribution of the 1,797 individual 
measurements within 742 groups with accurate distance measurements. Top right. The distribution of % 
errors in distances and peculiar velocities in the observed catalog. The median error of the CFl catalog is 
13% in both distance and peculiar velocity. Bottom left. Percent error in distance as a function of distance. 
Bottom right. The absolute errors on peculiar velocities increase with the depth of the catalog. 
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Fig. 2. — Top: the V8K redshift catalog containing 30,126 redshifts displays the structure in the distribution 
of galaxies within a box of radius 8000 km/s (black dots). The distribution of 1,797 galaxies with measured 
distances in Cosmicflows-1 is superposed in red. The Galactic plane lies at supergalactic coordinate SGY=0 
and accounts for a band deficient in sources. Bottom: an AitofF all sky projection of the 742 grouped peculiar 
velocities in Cosmicflows-1 with galactic extinction over-plotted in grey shading. 
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SGX (Mpc) SGX (Mpc) 

Fig. 3. — Reconstruction of the mock catalog, sampled from a BOX160 constrained simulation of the 
local universe. The upper row shows the density and velocity fields of the simulation, smoothed at 1 (left) 
and 10 h~^Mpc (right). The circle show the zone where a mock catalog of peculiar velocities, similar to 
cosmicflows-1, was created (the data-zone). The lower row shows the WF reconstruction within the data- 
zone (left) and at a larger scale than the data-zone (right). The contours represent the logarithm of the 
recovered density field normalized by the mean matter density of the universe (positive densities only), and 
the arrows represent the recovered velocity field. The length of the arrows are proportional to the amplitude 
of velocity. Note that the reconstruction is made only from the radial peculiar velocities of the mock data 
drawn from the central sphere of radius R = 30h~^Mpc (the data-zone), with a poisson distribution around 
15% error. There is excellent agreement within the 30h~^Mpc ring between the smoothed input (top right) 
and WF reconstruction at the same smoothing scale (bottom panels at two scales). 
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Fig. 4. — Aitoff all-sky maps in supergalactic coordinates of reconstructed densities and velocities in shells 
within 40 Mpc. From top to bottom: shells centered at h/h = 6.7 Mpc, IQ/h = 13.5 Mpc, Ib/h = 20.3 Mpc, 
and 25//i = 33.8 Mpc where h = 0.74. Left column: the color map and contours give the WF linear density 
field reconstructed from the Cosmicflows-1 catalog of distances. The black dots identify galaxies from the 
V8K catalog within a shell of i? ± 3.5 Mpc. Right column: the color map gives the reconstructed radial 
velocity field in the same distance shells. 
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Fig. 5. — Orthogonal views of WF reconstructions of local structure. Slices of ±10 Mpc. SGX-SGY slice is 
centered at SGZ=0 and SGY-SGZ slice is centered at SGX=0 while SGX-SGZ slice is displaced to SGY=16.9 
Mpc to coincide with the Virgo Cluster. At left, over density contours are spaced at 0.3 intervals on a linear 
scale. At right, colors represent over and under de:j^ities. Individual galaxies from the V8K catalog are 
plotted on the right. The local part of reconstructed 3D peculiar velocities are shown by black arrows on 
the left. Observed peculiar velocities are represented by red (outward) and blue (inward) arrows in the top 
two panels that include our location at the origin. 
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Fig. 6. — Aitoff all-sky maps in supergalactic coordinates of the radial part of the reconstructed divergent 
velocities in shells within 40 Mpc. From top to bottom: shells centered at 5/h = 6.7 Mpc through 30//i = 40 
Mpc where h — 0.74. The color map and contours give the WF divergent field reconstructed from the 
Cosmicflows-1 catalog of distances. The black dots identify galaxies from the V8K catalog within shells of 
i?±3.5 Mpc. 
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Fig. 7. — Comparison of the reconstruction with the actual galaxy distribution in the Local Universe at 
a larger scale. Top row: WF reconstructed density contours and full velocity vectors evaluated in a box of 
±80 Mpc sides. The two slices are ±10 Mpc thick with the top left panel centered at SGY=0 and the top 
right panel centered at SGY=16.9 Mpc, the SGY depth of the Virgo Cluster. The green dots are galaxies 
from the V8K redshift catalog with the slices. The over density contour spacing is 0.3 plotted in linear scale. 
The mean density contour is plotted as a red solid line. Bottom row: the reconstructed WF velocity field 
is decomposed into its local (left) and tidal (right) components where the local velocities are a response to 
the WF density reconstruction within 80 Mpc. The decomposed velocity field is represented by flow lines. 
Note that the local component shows the expected close correspondence with the density field, with flow 
converging to the high density peaks. 
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Fig. 8. — Separation of the WF velocity reconstruction into local and tidal components in two projections. 
At left, an SGX-SGY plane at SGZ=0; at right, an SGX-SGZ plane at SGY=16.9 Mpc (coincident with 
the structure containing the Virgo Cluster). On top, the full velocity reconstructions. These vectors are 
decomposed into (middle) the local components intMfced by mass within 80 Mpc and (bottom) the tidal 
component, the residual after subtraction of the local velocities from the full velocities. The concentric 
circles delineate the 40 Mpc extremities of the data-zone. 



